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Increasing experimental and theoretical evidence points towards a dual nature of the 5/ electrons 
in actinide-based strongly correlated metallic compounds, with some 5/ electrons being localized and 
others delocalized. In a recent paper (PRB xxx, 2004), we suggested the interplay of intra-atomic 
correlations as described by Hund's rules and a weakly anisotropic hopping (hybridization) as a pos- 
sible mechanism. The purpose of the present work is to provide a first step towards a microscopic 
description of partial localization in solids by analyzing how well various approximation schemes 
perform when applied to small clusters. It is found that many aspects of partial localization are 
described appropriately both by a variational wavefunction of Gutzwiller type and by a treatment 
which keeps only those interactions which are present in LDA-I-U calculations. In contrast, the en- 
ergies and phase diagram calculated within the Hartree Fock approximation show little resemblence 
with the exact results. Enhancement of hopping anisotropy by Hund's rule correlations are found 
in all approximations. 
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I. INTRODUCTION 



Many fundamental phenomena in solid-states physics 
can be considered as resulting from the competition of 
the repulsive Coulomb interaction trying to keep elec- 
trons as far from each others as possible and the kinetic 
energy term trying to delocalize electrons as much as pos- 
sible. One of the most prominent examples is the Mott- 
Hubbard transition from a narrow-band metal with delo- 
calized electrons to an insulator with localized electrons. 
Another one is band magnetism where the occupation 
of majority and minority bands differ dramatically. The 
resulting magnetic moment is typically aligned along cer- 
tain crystallographic axes by an anisotropy energy which 
is much weaker than typical electronic energies given 
by the bandwidths. The recently suggested scenario 
of a dual nature of the 5/ electrons in some uranium- 
based heavy-fermion compounds combines features of 
both these examples: There is experimentalii2i2i^i^*Sii 
and theoretical evidenceSiSiSiiSiiiiiSiiii^ that in these 
materials some of the 5/ electrons with j = 5/2 are local- 
ized and others are not. Which of the 5/ orbitals form lo- 
cal moments and which form bands will depend on the ge- 
ometry of the crystallographic environment of the partic- 
ular actinide site, which leads to small differences of the 
hopping (hybridization) matrix elements^ Band struc- 



ture calculations augmented by a mass-renormalization 
of the 'free' 5/ electrons by crystal field excitations of 
the 'localized' 5/ electrons lead to very good agreement 
of both areas and masses with dHvA data in UPts and 
UPd2Al3iSiii4iS In both cases, the best agreement was 
achieved treating the jz = ±3/2 orbital as itinerant and 
including two / electrons in the ionic core. Calculations 
within the so-called LDA-I-U schemei^ for XiVA^-xP^x 
also suggest two fully occupied / electrons per U atom, 
however with jz — -1-5/2 and jz — -f3/2 having occu- 
pations very close to unity^ A treatment within the 
self-interaction corrected local density approximation— 
(SIC-LDA) found for UPta that the configurations of 
lowest energies have a local moment iS The same au- 
thors found for UPd2Al3 a groundstate characterized by 
the coexistence of localized / ^ and delocalized Uranium 
/ electrons. Rare-earth metals and their sulfides, i.e., 
4/ electron systems, were also studied within the SIC- 
LDA and discussed in terms of localized and delocalized 
orbitals. However, the physical reasoning is quite dif- 
ferent from the present one. 

In a recent paper fi4 we suggested the interplay of the 
intra-atomic correlations described by Hund's rules and 
a weakly anisotropic hopping as a possible mechanism 
leading to a separation of the 5/ orbitals into band- 
forming and local-moment-forming ones. The underlying 



2 



idea can be described as follows. According to Hund's 
rules, atoms with three or two / electrons have angular 
momenta of J = 9/2 with 10-fold degeneracy and J = 4 
with 9-fold degeneracy, respectively. Consider now tran- 
sitions M{jz) — {^J'z\'^j A^J z) between any given pair 
of those degenerate states due to an annihilation oper- 
ator Cj^ stemming from the kinetic energy term of the 
Hamiltonian. Obviously, these matrix elements depend 
strongly on jz , leading to much larger transition ampli- 
tudes for some jz than for others. Turning this argument 
around, the following scenario evolves: The largest prod- 
uct of the matrix elements M{jz) multiplied by the over- 
lap tj_, between the corresponding atomic wavefunctions 
of neighboring sites determine which pairs of J = 9/2 
and J = 4 configurations are realized, thereby maxi- 
mizing the corresponding energy gain. In most cases, 
the other hopping processes are then surpressed. This 
scenario is supported by exact diagonalizations of small 
clusters presented in Refs. [141 andll9t 

This mechanism differs from other routes to partial 
localization and a dual nature in strongly correlated sys- 
tems with degenerate d or f orbitals which were sug- 
gested, e.g., for manganites. For the latter, the physics is 
different: The degeneracy is split by strong crystal field 
with the lower t2g orbitals forming a localized high-spin 
state while the eg electrons are delocalized. 

The numerical effort of exact cluster diagonalizations 
grows rapidly with system size. Any quantitative de- 
scription of bulk material will necessarily involve some 
approximations. Over the last decades, considerable ex- 
perience has been accumulated about the capability of 
various approximation schemes to describe many-body 
phenomena such as the above mentioned Mott-Hubbard 
transition and band magnetism^ The purpose of the 
present paper is to provide a first step towards a mi- 
croscopic description of partial localization in solids by 
analyzing how various approximation schemes perform 
when applied to the cluster model of Ref. [l4| The latter 
is briefiy summarized in Sect.^ In Sect. IIIII and IIVI two 
widely used approximation schemes are applied. In the 
first, the Coulomb interaction is approximated by its di- 
agonal part (direct term plus exchange term) in analogy 
with a modified local density approximation to density 
functional theory (LDA-I-U). The second scheme is mean- 
field theory, also known as the Hartree-Fock approxima- 
tion (HF). In Sect.0 we study a variational wavefunc- 
tion similar to that proposed by Gutzwillcr in his famous 
paper on the fcrromagnetism of transition metals. '^^ We 
end in Sect. IVll with a comparision of the various approx- 
imative treatments and with some tentative conclusions 
and speculations about the performance of these approx- 
imation schemes in bulk materials. 



II. THE MODEL HAMILTONIAN 

The simplest model describing the interplay between 
intra-atomic correlations of / electrons and weakly 



anisotropic hopping is considered in Ref. 13 

H = Hcoui + Ht + Hh. (1) 

The Coulomb interaction is written in terms of the usual 
fermionic operators (n) creating an electron on site n 
with angular momentum j = 5/2 and z component jz as 

-ffcoui = (2) 

IJ2 J2 Uj^,j.,3.si^. 4.1 ('^)42 ('^)s.4 (")s.3 (") ' 

with Clebsch Gordon coefficients C " and Coulomb pa- 
rameters Uj 



U. 



JziJz2 \U\J 



ETJ f^JJz f^JJz 



(3) 



5/2,],4 



The kinetic energy operator describes hopping (hy- 
bridization) between all pairs of neighboring sites (nrn) 



Ht = - tj^(^c+{n)cj^{m) + h.c. 



(4) 



A ferromagnetic or antiferromagnetic external or internal 
field 



Hh = - 



K jz c\ in) c ■ in) 



(5) 



is included as well. 

For the cases of two sites occupied by 5 electrons and 
for three sites and 8 electrons, iterative sparse-matrix di- 
agonalization easily yields the respective ground states 
and some low-lying eigenstates. The sizes of the Hamil- 
ton matrices are 600 x 600 and 18, 000 x 18, 000, respec- 
tively, provided that site configurations are limited to / ^ 
and /3 states, and 792 x 792 and 61, 758 x 61, 758 if this 
restriction is liftedii^ For the sake of a transparent pre- 
sentation, we focus in the present paper on the two-site 
model and stay in the 600 dimensional f^-f^ subspace. 

The Uj parameters are more robust to changes in 
the chemical composition or the application of pressure 
than the the tj^ . Thus, we use Uj parameters appro- 
priate for UPts, Uj=4 = 17.21el/, Uj=2 = 18.28el^, 
Uj=o = 21.00el^)^, and study the phase diagram in the 
tj^ space. As in Ref. IT^ we choose = ^5/2 in the fig- 
ures. The phase diagrams of the model Hamiltonian (Q- 
Q obtained by exact numerical diagonalization is repro- 
duced in Fig. n 



III. DIAGONAL INTERACTION MATRIX 
ELEMENTS 

Almost all electronic structure calculations for ex- 
tended systems and large molecules are performed ei- 
ther in the framework of the local density approximation 
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FIG. 1: Phase diagram for the f-f model of Ref.[T3 in the 
^1/2-^3/2 plane (^5/2=^1/2)) derived from the total magneti- 
zation (2 component of the total angular momentum) of 
the ground state for infinitesimal field hn =0^. Numbers 
indicate Jz, which is a good quantum number. 

(LDA) or using the Hartree-Fock approximation. The 
former accounts for correlations only in a rather crude 
way, while the latter neglects them altogether. Vari- 
ous recent papers suggested the inclusion of strong lo- 
cal correlations via an on-site energy-term of Hubbard 
type, ~ Uhjfik, in LDA calculations within the LDA-I-U 
schemei& or the related SIC-LDAiSiifliiS This is a promis- 
ing route for some groundstate properties of highly cor- 
related systems. 

The key approximations of the LDA+U scheme are on 
the one hand the LDA itself and on the other hand the 
inclusion of only the diagonal interaction matrix elements 
(direct and exchange) 

(6) 

The fact that this approximation corresponds to the 
LDA-I-U approach is best seen by inserting it into the 
Coulomb Hamiltonian |(2J) and rewriting the latter in 
terms of haj = jC^ j as 

^Coul = XI ^{^^j'z\U\jzj'z) naj.naj'^- (7) 
1=1,2 

Henceforth, we will refer to Eq. |(7J) as 'diagonal approx- 
imation.' 

Figure 121 presents the resulting phase diagrams in the 
plane of hopping parameters ti/2 = ^5/2 A 
visually striking difference to the result of the full cal- 
culation. Fig. ^ is the complete absence of the isotropic 
line ti/2 — ^3/2 — ^5/2 ill the phase diagram of the di- 
agonal approximation to the Hamiltonian. This reflects 
the fact that the condition JBJ of diagonality of U with 
respect to z-component indices breaks rotational invari- 
ance. The z component of the total angular momentum 
remains a good quantum number. We note that the diag- 
onal approximation overemphasizes the parameter space 
of phases with either strong ferromagnetic or strong an- 
tiferromagnetic alignment along z, i.e., Jfz = 15/2 and 
Jz = 1/2. The simplified treatment quite correctly 



FIG. 2: Phase diagram within the diagonal approximation 
Q in the ti/2-i3/2 plane {ti/2 ~ ^5/2), derived from the total 
magnetization jTz of the ground state for infinitesimal field 
hn = 0"*". In comparison with Fig. note the different plot 
ranges. 

predicts the ferromagnetic 'high-spin' state J'z — 11/2 
for ^ ti/2 = ^5/2- In general, the approximation 
fares worst in the case of small hopping anisotropy. The 
jTz = 5/2 phase is missing, while a J'z — 3/2 state is 
predicted only for rather strong hopping. 

IV. HARTREE-FOCK APPROXIMATION 

The Hartree-Fock approximation is one of the most 
widely used calculational approaches for interacting sys- 
tems. For its implementation, one uses the single-particle 
density matrix 

(8) 

written as sum of A'tot HF orbitals, which are superposi- 
tions of atomic orbitals described by the 

a=l,2 j^=-5/2 

We impose no restrictions on the single-particle states, 
which are arbitrary combinations of angular momenta. 
The system of HF equations isS^ 

3' 

with 

"2? = E ^^mfk') (11) 

kk' 

and a, a denoting the two different sites. The two-particle 
matrix element {jk\U\j'k') comprises both the direct 
Hartree and the Fock contributions. Note that it is not 
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restricted to two pairwise equal indices in the jz basis. 
The total energy is obtained from 



total 



2^^ 

jj' 



^1 



(HF) 
a,]]' 



(12) 



where double-counting of the interaction energy is 
avoided in the usual wayi2^ 

In general, the HF solutions do not have the total 
as a good quantum number. Therefore we cannot always 
use different jTz's to classify different phases. However, 
one still finds discontinuous jumps either of the expecta- 
tion value (J7z) or of its derivatives with respect to tj^ or 
h and we use these to define phase boundaries. Moreover 
for some phases with small kinetic energy, self-consistent 
solutions (^("^^ with good quantum numbers are found 
and, thus, Jz is a good quantum number in those cases. 
This applies to the phases with Jz = 15/2 and 17/2. For 
larger hopping parameters, one finds in general contin- 
uous variations of {Jz)- We denote phases as 'AM' or 
'FM' indicating antiferromagnetic or ferromagnetic spin 
correlations, respectively. Half-integer-valued subscripts 
give the approximate {Jz) expectation value found for 
tj^ parameters not too close to the phase boundaries. 



and the magnetic energy. For not too strong hopping, the 
two sites (anti-) align their angular momenta in a (anti- 
) ferromagnetic field configuration, see left (right) column 
of Fig. H 
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FIG. 3: Phase diagram in the ti/2-i3/2 plane (with = 
^5/2) in the Hartree-Fock approximation with ferromagnetic 
(right), vanishing (middle), and antiferromagnetic (left) fields. 
The phases are chararacterized from the z components of the 
angular momentum and the state occupation and labled ac- 
cording to the notation of Fig.0|and|^ 

Overall, the Hartree-Fock ^1/2-^3/2 phase diagram in 
Fig.Elshows only weak resemblance with the one from the 
exact calculation. The isotropic line is a phase boundary 
at zero field. 

Figure 0] illustrates the occupation pattern of four im- 
portant phases occuring in the phase diagrams shown in 
Figs. I1I3I a nd [7| of the present work as well as Fig. 5 
of Ref. I14I A large angular momentum as implied by 
Hund's rules is obtained if electrons occupy preferren- 
tially states with either large positive or large negative 
angular momentum projections jz - of course subject to 
the Pauli principle. Occupation of states with large posi- 
tive or negative angular momentum projections is favored 
also by local (molecular) magnetic fields. In particu- 



lar, the antisymmetrized product states 

.5 _i_3 



1-1-2: 



± |, minimize simultaneaously the Coulomb energy 



FIG. 4: Configurations relevant for weakly anisotropic hop- 
ping (hybridization) in the presence of ferromagnetic (FM, 
left column) and antiferromagnetic (AM, right column) fields 
for tz/2 < ti/2 = ^5/2 (top row) and ti/2 = < h/2 (bot- 
tom row). Rectangles symbolize the orbitals on both sites, 
with black areas presenting their occupation. The thickness 
of horizontal arrows between orbitals symbolizes their contri- 
butions to the kinetic energy gain. The arrow-shaped back- 
ground emphasizes the role of the fields. Note the charge 
disproportionation < An/ < 1 in the antiferromagnetically 
aligned states. 



Two crucial differences between FM and AF fields are 
obvious: (i) As drawn, in the FM case both sites have an 
/ occupation n/ = 2.5, whereas in the AM case a small 
charge disproportionation An/ < 1 can be present. The 
AM configurations of Fig.^are, of course, degenerate due 
to the symmetry with respect to a simultaneous exchange 
of site index and jz <-> —jz ■ 

[a) The FM and AM cases differ with respect to 
a possible quantum phase transition when going from 
^3/2 < ii/2 to ^3/2 > ti/2 (compare in Fig. 2]thc upper 
row with the lower row). In the AF case, with increasing 
ti/2 the jz = 1/2 electron starts to fluctuate more be- 
tween the sites, whereas the jz = 3/2 electron fluctuates 
less. This is a continuous cross-over, which however takes 
place rather suddenly because jz — 1/2 fluctuations im- 
ply charge fluctuations which then suppress jz =3/2 fluc- 
tuations, and vice versa. In a level diagram such as the 
inset of Fig. 4 in Ref. !l5, this corresponds to an avoided 
crossing. 

In contrast, in the FM case one of jz = 1/2 or jz = 3/2 
is not fully occupied. The transition is not a cross-over, 
but a simple level crossing with discontinuous change of 
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Jz and orbital occupation. This difference between sharp 
discontinuities on the FM side and smooth cross-overs on 
the AF side wiU be seen repeatedly in this work. We will 
see below that within the Hartree-Fock approximation 
the states FM15/2, FMyjji, and AM^ji dominate indeed 
in the phase diagram for finite fields and not too large 
hopping. 

One further cautionary remark is in place: The rep- 
resentation in Fig. ^ shows only the occupation pattern, 
but does not reflect the correlations between different 
configurations. 

Strongly anisotropic hopping i3/2 ^ ^1/2 lead in the 
full solution to FM15/2, i-c. to ferromagnetic correlations 
and single occupancy of the binding j^. +§ state. The 
corresponding wavefunction 



^3/2 



(a) 



^3/2 



(&) 



^5/2 



(6)c+ 2(6)10) 
(13) 

is a simple Slater determinant. Thus, it is not surprising 
that HF works well for this case. The parameter space 
of the FM15/2 phase is, however, much larger in the HF 
approximation than in the exact calculation. It shrinks 
with increasing antiferromagnetic field and increases with 
increasing ferromagnetic field. 

Outside the ferromagnetic phases in Fig. (SJ the corre- 
lations are antiferromagnetic and (jTj) < 3/2. A kink 
separates a phase with (jTz) ~ 1/2, visualized in Fig. ^ 
as AMi/2, from a phase with (J7z) ~ 1 at or near the 
isotropic line and (jTz) ~ 3/2 at larger values of ^3/2 
which is denoted by AM3/2 in Fig. [3 In a small bubble 
bordering the ferromagnetic FM]^5/2 phase, a phase exists 
which we denote by AM3y2, see Fig.|Sl It can be consid- 
ered as the symmetrized version of AM3/2- The bind- 
ing linear combinations with jz = +3/2 and jz = —3/2 
are both half occupied and contribute equally to the ki- 
netic energy gain by hopping (hybridization). Strongly 
anisotropic hopping with ti/2 ^ ^3/2 leads to AM1/2, i-e. 
to an antiferromagnetic arrangement with broken left- 
right symmetry, i.e., < An/ ^ 1. 

A region of maximal ferromagnetic moment (FM17/2) 
is found for a ferromagnetic field which is still weak com- 
pared to the bare hopping strengths. For larger fields, 
the ferromagnetic phase FM5/2 with 'intermediate' spin 
is found close to isotropic line. In this phase, all five 
electrons contribute to the kinetic energy gain. 

Figure [SI shows also a phase denoted by FMii/2j which 
is of overall ferromagnetic nature but dominated by ^3/2 
hopping. It appears in Fig. and |2] as well as in Fig. |S| 
below. 



A. Partial localization 

In order to quantify the degree of localization of a given 
jz orbital by local correlations, we study the ratio of the 
jz-projected kinetic energy Tj^ and the bare matrix ele- 
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FIG. 5: Configurations relevant for large hopping: The upper 
left panel (FM5/2) shows an overall ferromagnetic arrange- 
ment with all electrons contributing to the kinetic energy gain. 
The upper right panel (AM3/2) presents an overall antiferro- 
magnetic arrangement dominated by £3/2 hopping which is 
symmetric with respect to the site index. 1^/2 hopping domi- 
nates also in the lower row. In the lower left panel (FM11/2), 
the overall arrangement is of ferromagnetic type, whereas in 
the lower right panel (AM3/2) it is of antiferromagnetic char- 
acter. As in Fig. |1] black areas symbolize the orbital oc- 
cupation and the linewidth of the arrows symbolize the cor- 
responding contributions to the kinetic energy gain. Field 
directions are not indicated, in order to emphasize the role of 
the hopping for the formation of these states. 



ment t. 
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E (*g«l(4,.Wc±,.(m) + /i.c.)|vE'g,), (14) 



with —2 < Tjjtj^ < +2 for two sites. A small ratio 
Tj^ /tj^ indicates partial suppression of hopping for elec- 
trons in the ztjz orbitals. For all cases considered in 
Ref. 113. we found small ratios of Tj^ / tj^ for the smaller 
tj^ , and values close to one or two for the dominant hop- 
ping matrix element. 

In Fig. 1^1 we study the same quantity for the Hartree- 
Fock solutions. Again if ^ ^1/21^5/21 ratios close to 
one are found for 13/2/^3/2, but much smaller values for 
(T'1/2 + 1^1/2, and vice versa. At first glance, this 
may look surprising since the Hartree-Fock method can 
not describe the correlations underlying the dual nature 
of the / electron in our model system. A closer look 
reveals that the correlations are simulated by symme- 
try breaking, i.e. by driving the corresponding occupa- 
tions close to zero or two. In contrast, the full solution 
shows a reduced hopping also for intermediately occu- 
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pied orbitals. We conclude that HF mimics correlations 
by overemphasizing occupation differences. Thus, fur- 
ther studies are needed to gauge the usefulness of HF for 
bulk systems with partially localized 5/ electrons. 




0.0 0.5 1.0 1.5 0.0 0.5 1.0 1.5 
ti/2 = t5/2 (eV) ti/2 = t5/2 (eV) 



FIG. 6: Ratio (a) T^/2/tz/2 and (b) (Tj/a + T5/2)Ai/2 in 
the Hartree-Fock approximation in the 1^/2-^3/2 plane (with 
ti/2 = ^5/2 and hn ~ 0^) given as gray-scale plot (black: 0, 
white: 1.5). Phase boundaries of the middle panel of Fig. |21 
are included for orientation. 



As an aside, we remark that numerical convergence in 
solving the restricted equation H15|l is obtained in a small 
fraction of the time needed for solving the HF equations 



h " = 0.01 eV h = h = 0.01 eV 




1 2 30 1 2 30 1 2 



'1/2=^5/2 (eV) t„2 = t5,2 (eV) ti,2 = t5,2 (eV) 



FIG. 7: Hartree-Fock approximation with variational space 
restricted to eigenstates (UHF): Phase diagram in the - 
i'i/2 plane (fi/2 = ^5/2) with ferromagnetic (right), vanishing 
(middle), and antiferromagnetic (left) fields. The phases are 
chararacterized from the z components of the angular mo- 
mentum and the state occupation, cf. Fig. 01 and |S] 



B. Hartree-Fock with eigenstates (UHF) 

If the variational HF space is reduced further by re- 
quiring j^-diagonal density matrices and self-consistent 
potentials 

^a,r^ - E^-^^fc 0-fc|V^b-fc>, (15) 

fe 

the z axis plays a special role (quantization axes), as 
it did in Sect. IIIII and the spherical symmetry is bro- 
ken. Following the usual nomenclature, we refer to this 
approximation which allows different orbitals for differ- 
ent angular momenta but prescribes angular momenta 
as orbital quantum number, as the unrestricted Hartree- 
Fock (UHF)2^ - even though the wavefunction space is 
restricted in comparison to Eq. (j^. The resulting phase 
diagrams. Fig. [T] are not identical, but close to those 
of the full HF solution in Fig. |31 Not surprisingly, the 
isotropic line does no longer play a special role. The 
corresponding phase boundary has moved to higher tj,/2 
values, a consequence of the variational restriction fa- 
voring the maximally antiferromagnetic alignment with 
Jz = 1/2. 

Note that Jz is now a good quantum number. In par- 
ticular, the border between AM3/2 and AM1/2 is here 
easily identifiable, the kink being replaced by a jump. A 
noticeable difference between the HF and the UHF phase 
diagram is that AM^^^^ i.e. the symmetrized version of 
can not be found within the UHF approxima- 
tion, because all six values are partially occupied in 
this phase, cf. Fig. [S] 



C. Charge disproportionation 

The Hartree-Fock phase diagram is surprisingly rich 
with several phases showing charge disproportionation, 
i.e. a finite expectation value < An/ < 1 of the dif- 
ference of site occupations. Note in particular the phase 
boundary in Fig. 13 between the closely related phases 
AM3/2 and AM3^2 at vanishing or staggered field. 

The exact groundstates of the small cluster does not 
show charge disproportionation, because a suitable su- 
perposition of the two symmetry-related configurations 
will lead to a finite energy gain, completely analogous to 
the absence of ferromagnetism in small clusters at van- 
ishing external fields. However, in both cases the en- 
ergy gain can be extremely small if the configurations 
are sufficiently different. This feature of the numeri- 
cally exact small cluster solution will often not survive 
the thermodynamic limit. In contrast, mean-field ap- 
proaches easily lead to symmetry-broken solutions even 
for small clusters. This is so, because they reduce charge 
fluctuations and therefore simulate the effect of electron 
correlations. When symmetry breaking is experimen- 
tally observed in bulk samples it is overemphasized in 
a Hartree-Fock treatment. Thus, we take the trend of 
the HF solution towards charge disproportionation only 
as a hint that it might be worthwhile to look for, e.g., 
phase transitions driven by charge ordering in systems 
with partially localized 5/ orbitals. 

V. GUTZWILLER'S WAVEFUNCTION 

The exact groundstates found in Ref. ^3 often have a 
large overlap with a Gutzwiller-type wavefunction^i of 
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the form 



(a) 



h=0 



h=0.01eV 



(16) 



where the projector P/a/a retains only local and 
configuration. The single-particle orbitals generated by 
the are subject only to the orthogonality requirement. 
In most actual applications of the Gutzwiller ansatz, 
simple trial wavefunction are used for the single-particle 
states, e.g., the non-interacting, possibly polarized, Fermi 
sea. Here, we use the ansatz that each orbital has jz as 
a 'good quantum number,' i.e. 

5 

\^a)=Vpp n (cos0,c+„)(a) + sin0,c+„,(6))|O), 

(17) 

with variational parameters . Both, the numerator and 
the denominator of the energy expectation value 



E[{ji-'^},{9^}] 



(*G I *g) 



(18) 



Note that 



are polynomial expression in cos 9^ and sin 9,j 
in the case of double occupancy jz"^^ = ji 
sponding variational parameters 9,^ drop out of Eq. (|18|l . 

We performed numerical minimizations with respect to 
both the occupation pattern {ji'^^ } and the set of param- 
eters {6^}. In practice, the expectation value l(TH)l is ob- 
tained explicitly from the Hamilton matrix in the many- 
particle site-spin representation (600x600) in terms of 
9^. Then Mathematica's implementation of Brent's 
algorithm is used for 1-dimensional, 3-dim., and 5-dim. 
numerical minimization in the case of 2, 1 and doubly 
occupied angular momenta jz, respectively. Even the 
five-dimensional minimum search is fast. Unfortunately, 
many occupation pattern have to be tried and different 
random initial conditions yield different answers. This 
makes an additional loop over various initial conditions 
necessary. 

The obtained phase diagram is reproduced in Fig. |SJi. 
Again, J'z is a good quantum number because the cor- 
responding operator commutes with Vpp. Both panels 
show a transition from a strongly polarized J7z = 15/2 
phase to the state of lowest possible polarization, i.e., 
Jz = 1/2. A smah pocket with Jz = 17/2 is found at fi- 
nite field. The phase diagram for small antiferromagnctic 
fields (not shown) is close to that for vanishing field with 
the high-spin phase Jz = 15/2 being slightly reduced. 
As in the exact solution and in all other approximation 
schemes, the energy varies continuously (not shown). 

The phase diagram of the Gutwiller wavefunction 
is remarkably close to the exact diagonalization of 
the Hamiltonian keeping only diagonal matrix elements 
('LDA-hU') in Sect.llTIl We wiU see below that this is true 
for the total energies as well. This should not come as 
an surprise, because both approaches effectively reduce 
charge fluctuations but do not include angular correla- 
tions equally well. 
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FIG. 8; Phase diagram for the Gutzwiller wavefunction in 
the ti/2-t3/2 plane (with ti/2 = ^5/2) for vanishing (left) and 
finite (right) ferromagnetic field, (a) Phases are character- 
ized by their expectation value, (b) Ratio 73/2/^3/2 as 
gray-scale plots (black: 0, white: 1.4) and (c) the sum of the 
corresponding ratios for = 1/2 and 5/2 (black: 0, white: 
2.4). In comparison with Fig.|5| note the different plot ranges. 



A. Partial localization 

In Fig. [SId and IS]:, we study the quantity Tj^/tj^ for 
the Gutzwiller type ansatz H17|l . If either t3/2 ^ ^1/2 or 
^1/2 S> ^3/2, we find again a clear suppression of the non- 
dominant hopping. It is interesting that in the Jz = 1/2 
phase one can go from one case to the other without 
crossing a phase boundary. This is interesting in so far as 
a Gutzwiller wavefunction describes this important point 
without having to introduce phase transitions by symme- 
try breaking. 



VI. SUMMARY AND OUTLOOK 

We have studied the interplay and competition of intra- 
atomic correlations and an anisotropic kinetic energy op- 
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FIG. 9: Total energies for various approximation as function 
of ~ t5/2 at fixed t3/2 = 1.65 eV: (exact) numerical di- 
agonalization of the 600 x 600 matrix, see Fig. (diagonal) 
diagonalization keeping only the diagonal interaction matrix 
elements, see Fig. |21 (HF) Hartree Fock approximation, see 
Fig- 0] (UHF) Unrestricted Hartree Fock approximation, see 
Fig. |7| (Gutzwiller) variational wavefunction 1171 , see Fig. |H| 

erator by applying three different approximative many- 
body treatments to the model of Ref. Q. None of them 
leads to fully satisfactory results. But, the diagonal ap- 
proximation inspired by the LDA+U approach, i.e. keep- 
ing only the matrix elements gives a reasonable good 
phase diagram with most important phases of the exact 
solutions being present. This finding should encourage 
further applications of LDA-t-U and related approaches 
such as SIC-LDA for groundstate properties of actinide 
heavy-fermion materials. We mentioned already in the 
introduction that some LDA+U calculations indeed show 
partial localization. However, one can not expect to 
reproduce within these schemes the small energy scale 
responsible for the heavy-fermion character of the low- 
energy excitations. 

The numerically easier UHF approximation with j^- 
diagonal orbitals agrees quite well with the result of the 
most general HF ansatz. However, both yield phase di- 
agrams with only weak resemblance with the one from 
the exact calculation. The substantial suppression of the 
subdominant hopping found in Fig. might be somehow 
fortituous. 

These observations are confirmed by the total energies 
given in Fig.|51for all approximation schemes along a line 



^3/2 = 1.65 eV in the ti/2-H/2 plane. The particular 
^3/2 value was chosen such that a large number of phase 
boundaries is crossed. Clearly, at some phase boundaries 
kinks of the total energy are seen (related to level cross- 
ings) whereas at others higher derivatives change more 
or less discountinuously (avoided crossings). 

We note again that for the exact numerical diagonaliza- 
tion and for the calculation keeping only state-diagonal 
interaction matrix elements we surpressed configurations 
/" with 71 7^ 2, 3, i.e., we worked in the 600 dimensional 
space rather than the 792 dimensional space. Giving 
up these restrictions lowers the corresponding energies 
slightly (not shown). 

All curves, except that of the exact solution start with 
the same value which is constant for a certain range of 
ti/2 and values. The corresponding wavefunction 
is the Slater determinant 1 4") 15/2, Eq- Ijl^l) . which com- 
pletely surpresses ti/2 and i5/2 hopping. In contrast, the 
exact ground state near = 1.65 eV, ti/2 — tc^/2 — 
OeV is not |^')i5/2, but has intermediate angular mo- 
mentum Jz = 5/2 and gains energy from all hopping 
processes. This explains the lower energy and its depen- 
dence on all hopping matrix elements. 

The solutions of the most general Hartree- Fock scheme 
(HF) and the case where the single-particle orbitals were 
restricted to be diagonal with respect to jz (UHF) differ 
only in small parts of the parameter space. 

While the present model was developed with actinide- 
based materials in mind, it seems worthwhile to examine 
the role of intra-atomic correlations in other materials, 
most noticeable in transition metal compounds. 

In conclusion, we find the solutions for the modified 
Hamiltonian containing the diagonal matrix elements 
only and the expectation value of the original Hamilto- 
nian with generalized Gutzwiller variational wave func- 
tion H17|l are quite close, because both approaches ef- 
fectively reduce charge fluctuations but do not include 
angular correlations equally well. The subtle angular 
correlations which determine the magnetic character are 
accounted for only in limiting cases. 
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